Forcing-dependent stability of steady turbulent states 
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We study the influence of tfie forcing on the steady turbulent states of a von Karman swirling 
flow, at constant impeller speed, or at constant torque. We find that the different forcing conditions 
change the nature of the stability of the steady states and reveal dynamical regimes that bear 
similarities with low- dimensional systems. We suggest that this forcing dependence may be an out- 
of-equilibrium analogue of the ensemble inequivalence^ valid for long-range interacting statistical 
systems, and that it may be applicable to other turbulent systems. 
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Introduction An intriguing property of statistical sys- 
tems with long range interactions is the ensemble inequiv- 
alence: a solution in the microcanonical (constrained) 
ensemble is not necessarily a solution in the canonical 
(unconstrained) ensemble [l|, This property traces 
back to the non-additivity of energy, and is reflected by 
pathological behaviours of the entropy (that can be non- 
concave) or the heat capacity (that can become negative). 
Ensemble inequivalence has been observed and studied 
in a variety of systems such as for example 2D Euler 
equations, [sl, Jj, Blume-Emery- Griffiths model [5], ran- 
dom graphs \&\. More recently, it has also been studied 
in a model describing the one-dimensional motion of N 
rotators coupled through a mean-field interaction, and 
subjected to the perturbation of an external magnetic 
field 0; it shows that this concept may also hold out of 
equilibrium, broadening its range of applicability. Nev- 
ertheless, a generalisation of such concepts to all out-of- 
equilibrium systems remains — at present time — an 
open question. We propose an extension of such tools 
to turbulent systems. They naturally display long-range 
interactions and are, by definition, far from equilibrium. 
Our experimental system is a "von Karman" experiment, 
in which a cylinder of fluid is stirred by two impellers 
with blades, producing fully-developed turbulence in a 
relatively small experimental device. Several features of 
equilibrium systems have already been described in this 
model experiment, such as steady states, predicted by 
the resolution of the axisymmetric Euler equations [8|, 
hysteresis [9], and spontaneous symmetry-breaking with 
diverging susceptibility ^lOj . 

In this letter, we will examine the stability of such 
steady states under two different forcing conditions, ei- 
ther imposing the speed or the torque to our impellers. 
The subject has attracted little attention, only focusing 
on the difference of power fluctuations JjJ, ll2| under both 
conditions. However, we can consider our forcings to be 
conjugate, as the product torque x speed controls the 
energy injection rate in our experiment. Switching from 
speed to torque control might then be seen as an analogue 



of switching from canonical to microcanonical ensemble; 
it is shown to alter the stability of the steady states previ- 
ously observed in [ol and to reveal interesting dynamical 
regimes. 

Experimental setup The von Karman flow is created 
in a Polycarbonate cylinder of radius R = 100 mm filled 
with water. The fluid is stirred by two 8-bladed impellers 
of radius 0.925i? separated from each other by a distance 
(blade-tip to blade-tip) 2h = lAR (see fig. [I]). Two in- 
dependent 1.8 kW brush- less motors can rotate the im- 
pellers by imposing torque - or speed - to the system. 
The fluid is confined in the cylinder by two balanced me- 
chanical seals under a P = 2.8 bar pressure to ensure 
minimum friction. Two torque sensors measure the ro- 
tation angle - and therefore the speed - and the torque 
applied on the two impellers, respectively (/i,Ci) and 
(/2,C2). Temperature control is enforced by an external 
water flow in two refrigeration coils installed behind each 
impeller. We define the Reynolds number of our exper- 
iment as Re = Illli^-h}^^ varying from 2 10^ to 5 10^ 
depending on the speed of the impellers. Furthermore, 
we define two main parameters, each corresponding to 
the control parameter of one type of experiment: 

• = the reduced impeller speed difference 

for the experiments where 2/o = /i + /2 = 8 Hz is 



imposed; 



C1—C2 



the reduced shaft torque difference for 
Ci + C2 = 3.45 Nm 



^ - C1+C2 

the experiments where 2Co 
is imposed. 



/ and C are conjugate quantities: for a given exper- 
iment, we either impose a torque C - respectively an 
impeller speed / - to measure the impeller frequency / 
- respectively the torque C. The product of these two 
quantities also defines the amount of power injected in 
the experiment, V = 27r fC. 

Speed imposed For speed- imposed experiments, a hys- 
teresis cycle which opens for Re > 4000 has been ob- 
served [9]. Fig. [2] depicts the torque-speed asymmetry 
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FIG. 1: Sketch of the 'VK2' experiment. Left: Impeher used 
for the experiment and blade profile: curvature a = 72° , blade 
height 0.2R. Right: section of the experiment with the T^tt 
and (Oz) symmetries. The refrigeration coils are not shown. 



relation 7(^)1 starting both impellers at the same time 
for very close to 0, the system remains steady in a 
'symmetric' branch called (s). A slight variation of on 
this branch triggers a global bifurcation with a dramatic 
increase of the mean torque ^{Ci + C2). Two new, 'bi- 
furcated' branches of the (7,6^) plane are created; (61) 
(mainly for ^ > and 7 > 0) and (62) (mainly for ^ < 
and 7 < 0) are exchanged by TZ^^ rotation. 

These states have been characterized both by stereo- 
scopic Particle Image Velocimetry (sPIV) and Laser 
Doppler Velocimetry (LDV): in (s), each impeller gen- 
erates one flow cell by Ekman recirculation. The two 
cells are separated by a shear layer, similarly as in (ol, [iq| 
for ^ = 0. In contrast, both bifurcated branches display 
only one circulation cell covering the entire system. PIV 
measurements have confirmed that (61) and (62) niean 
velocity fields are effectively images of each other by the 
IZ-j^ rotation. 

Once on (bi) and (62) branches, it is never possible 
to reach again the {s) state, which is thus only observed 
while starting the impellers at the exact same time. In 
addition, the (61) and (62) branches are hysteretic, (61) 
being for example stable for down to —0.15. Further- 
more, a notable 7 range [-0.089; -0.018] U [0.00; 0.055] is 
not accessible due to the separation of our three branches 

Torque imposed In contrast, imposing torque allows 
us to choose 7, assuming friction is negligible. For 
the values of 7 accessible by imposing speed, the states 
reached are steady and we clearly see that the mean ob- 
served 6 is consistent with speed-imposed experiments. 






FIG. 2: Torque asymmetry parameter 7 = c\+cl ^ func- 
tion of the speed asymmetry parameter = j^^j^- Speed 
experiments ( ) impose 9 to determine 7, whereas torque ex- 
periments ( ) do the opposite. While a hysteresis cycle is 
observed in speed experiments, new 'forbidden' branches are 
observed. These branches are multistable, mixing slow bifur- 
cated — (^1,2) — states (v), intermediate — (^1,2) — (^), and 
fast symmetrical state (s) (a). The hatched region represents 
the inaccessible 7 values imposing speed. 



We will then assume that they are identical. This result 
has also been confirmed by PIV measures, with torque- 
driven and speed-driven velocity fields almost identical 
when properly normalised in both {s) and {b) branches. 
However, 7 values from the forbidden region are now ac- 
cessible. In such regions, the system looses its steadi- 
ness and transits between multiple turbulent states. The 
multi stability is characterized by the emergence of mul- 
tiple peaks in the P.D.F. (probability density function) 
of Ofit)^ the 1.5 Hz low-pass filtered signal of 0(t). This 
filtering is required considering our speed measurements 
generate discrete values of the impeller speed /, and yield 
a robust density function when the filter cut-off frequency 
is changed. Therefore, three types of "attractor" states 
have been defined: (5), the high-speed state, with charac- 
teristics similar to (5), (61) and (62) which are low-speed 
states similar to (61) and (62), and two new {ii) and (22) 
intermediate states. All of these states may hold longer 
(> 10 sec.) than any characteristic time of our experi- 
ment for specific 7 values inside the formerly forbidden 
zone. These new states can be seen in fig. [2l while (5), 
(^1) (^2) states extend well beyond their speed-imposed 
counterparts, the {ii) and (^2) branches are new and can- 
not be observed otherwise. 

Decreasing 7 from a perfectly symmetric (s) (6 = 0) 
state, we can observe its influence on the temporal sig- 
nals, (cf. fig. [3]). First, the speed of the impellers re- 
mains steady. Then (fig. [31^6)), when 7 < 7ci ^ —0.049 - 
corresponding to the extremum of small excursions 
towards (^2) occur. For larger 7, the oscillations be- 
come Jarger and allow sometimes short jumps towards 
the (62) state (fig. [3fc)). In the region where all (s), 
(^2), (^2) states are observed (fig. Wid))^ some "steady" 
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FIG. 4: Shifted temporal signals of 60 randomly- chosen tran- 
sitions of a 2-hour experiment with 7 = —0.0891. We compute 
Ti by finding the minimum of \dtf2[ where 7 means the . signal 
has been numerically filtered at 1 Hz. (a), (c), respectively 
/2 and /i profiles for down transitions. (6), (d), respectively 
/2 and /i for t transitions. The thick white line represents 
in each sub-plot the rotation frequency averaged on all 195 
events of the experiment, 



FIG. 3: Temporal series of the impeller speeds fi ( ) and 
/2 (red) for various 7. (a), symmetric high-speed state (s) 
observed at 7 = —0.0164; (6), threshold of the irregular os- 
cillations (^2) with very small events for 7 = —0.0460; (c), 
(22) irregular oscillations for 7 = —0.0668; (d), multi stable 
regime showing (s), (22) and (62) events at 7 = —0.0891; (e), 
a single fast rare event in a quasi-steady slow (62) regime at 
7 = -0.0912; (/), slow (62) regime for 7 = -0.1049. 



(12) events, exceeding 10 sec (70 impeller cycles), have 
been observed. Decreasing 7 further affects this region 
by favouring (b) at the expense of {i) and (s). Therefore, 
for low 7 > 7c2 ^ —0.0920 (fig. [3l^e)), only rare events 
can drive the system to the (^2) state, and nearly no time 
is spent in the (s) state. Eventually, the system joins the 
(62) steady branch (fig.[l/)) for 7 < -0.099. 

Valuable information on our system dynamics can be 
found studying the temporal signals of the global quan- 
tities near the transitions [13]. We have therefore super- 
posed in fig. m the speed signals close to the transitions 
observed in fig. ^c): (s) (^1,2) is called down, and 
(^1,2) (s) up. Once the transition times accurately 
determined, a fair collapse of /i,2([Ti — 5 sec; + 3 sec.]) 
has been observed. It is therefore possible to deduce a 
general transition scenario for up and down transitions. 
The up transition seems to result from a correlated pos- 
itive fluctuation of the two impellers, whereas the down 
one suggests a degenerate (^2) oscillation in which a sud- 
den impeller correlation is triggered very close to r^. 

Eventually, the distribution of (/i,/2) are studied to 



characterise the attr actor states that emerge from figs. [2] 
and [3l In fig. [Sb), for 7 too small, only one peak ap- 
pears, which confirms the steady nature of the system in 
{s). For small asymmetries, small excursions that escape 
the attractor, corresponding to the previously described 
small (z2) oscillations can be found, exhibiting a new peak 
strongly deviating from the diagonal /i = /2- For higher 
asymmetries, the system fills a large part of the (/i,/2) 
plane, with three main maxima: (s) close to the diago- 
nal at higher (/i, /2), (62) for the off-diagonal low (/i, /2) 
regime. The third — (^2) — attractor is harder to see be- 
cause it is partially hidden by neighbouring zones repeat- 
edly crossed by unsteady events. It is located roughly 
between the right tip of the joint histogram and the (s) 
maximum. With this representation, one clearly observes 
that the mean system path is different for the down and 
up transitions: while the down transition seems to start 
"looping" in the vicinity of (s) before abruptly transit- 
ing to (62), the up transition reaches the right edge of 
the joint-PDF (/i > /2), near (^2) before joining the (s) 
state. 

The repartition of the height of the maxima in fig.[5](c), 
(d), (e) is obviously driven by 7, from almost-fully (?), 
(^2) to nearly-pure (62) with rare large transitions to the 
faster states. For nearly-pure (s), we clearly see (fig EJc)) 
that a large amount of excursions occur, contrasting with 
the nearly-pure (62) state (fig El^e)), for which no excur- 
sion is observed (though a few transitions are). In addi- 
tion, the position of these maxima is not fixed when 7 is 
changed: indeed, we recall for example that fig. [2] shows a 
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FIG. 5: Joint-probability density maps of the (/i,/2) val- 
ues (log-scale), based on the temporal signals of fig. [S] (a), 
steady symmetric high-speed state: 7 = —0.0164; (6), begin- 
ning of irregularly oscillating states towards an intermediate 
state, 7 = —0.0460; (c), typical irregular oscillations with 
large slowdowns: 7 = —0.0668; (d), multistability with three 
'most visited' states: 7 = —0.0891, ( — ) and ( — ) represent 
respectively the fig. |4] mean profile for down and up transi- 
tions; (e), rare events for 7 = —0.0912; (/), steady slow state: 
7 = —0.1042. The black dashed-dotted line represents the 
^ = condition. Both fast and slow states can be observed 
at ^ = considering the shape of the fig. [2] curve. 



clear Q variation of the steady (s), (61,2) maxima when 7 
is increased. Interestingly, while the (6) and (?) Q values 
increase with 7, the (i) position decreases with 7. 

Discussion Using global torque and speed measure- 
ments, we have characterised the response of the von 
Karman experiment to different energy injection mech- 
anisms. The two responses coincide in the range of pa- 
rameters accessible to the constant velocity forcing mode, 
corresponding to the hysteresis cycle previously reported 
by [qI. However, imposing the torque 7 in the zone which 
could not be reached in speed experiments generate new 
continuous "mean" branches connecting symmetric {&) 
and bifurcated (b) branches. However , these "mean 
states" are not representative of the dynamics observed: 
the P.D.F of may exhibit up to three peaks, corre- 
sponding to "intermittent" states. Two of them can be 
defined by continuity with the speed-imposed states: (s) 
and (b). The third state, (z), is never observed in any 
speed-imposed experiment. The study of the impellers 
velocity Vi{t)^V2{t) signals show typical excursions and 
transitions between our three states, similarly to [l6l.[l7|. 
while preliminary results on the distribution of (6) resi- 
dence time favour exponential Kramers-like escape times. 
It might then be possible to describe the evolution of 
our global quantities with a very siniple model of ground 
states with noise, as performed by [18|. Work is also in 



progress to characterize the geometry of the correspond- 
ing steady states through PIV measurements. 

Our results address two questions. First, from the 
point of view of statistical physics, the von Karman ex- 
periment allows a quantitative analysis of the influence 
of the energy injection on the response of an out-of- 
equilibrium system. Our forcings are conjugated with re- 
spect to the energy injection rate, equivalent in a steady 
case to en energy dissipation rate, or entropy production 
rate, and do not provide the same results: more "working 
states" are found when imposing torque. This is remi- 
niscent of what happens in systems with long-range in- 
teractions, between the microcanonical (constant energy) 
case and the canonical (constant temperature) case. For 
a particular set of parameters, the microcanonical ensem- 
ble exhibits a negative specific heat region = < Q, 
despite continuous variation of T. Canonical ensemble 
cannot sustain such Cy < otherwise energy fluctua- 
tions in the system generate destabilizing heat transfer 
from the thermal reservoir. Therefore, a brutal first-order 
transition is observed, similar to our bifurcations in speed 
control, whereas torque control allows a continuous tran- 
sition. More specific features point towards ensemble in- 
equivalence in our system: a new quasi-steady solution, 
(i), arises only in torque control. In addition, positive 
X regions are found to be identical, "ensemble equiva- 
lent", in both controls, whereas only torque control al- 
lows multiple portions of negative, "ensemble inequiva- 
lent", X~^{l) = 1^5 which cannot be found in speed con- 
trol. These results strongly recall those obtained with 
out-of-equilibrium rotators [Tj. 

Second, from the point of view of turbulence, it sets the 
problem of universality of the steady states, that appears 
to be rather sensitive to the energy injection mechanisms, 
at variance with traditional view of turbulence. The phe- 
nomenon we explore in the present letter could actually 
be present in other turbulent experimental systems: for 
example, turbulent Plane-Couette flows could be forced 
either with constant global stress (motor torque C) or 
global strain (speed /), Poiseuille flows have be studied 
either inxQosing a pressure difference AP or a mass flow- 
rate Qm [201 , and even Rayleigh-Benard experiments have 
been conducted under temperature- imposed (T) and the 
heat-fiux imposed {dtQ) conditions [2l|. While the two 
first examples also exhibit a "duality" of the forcing types 
{T^ = Qm^P = Cf), Rayleigh-Benard experiments do 
not respect this condition {V = dtQ). It is difficult to 
understand what would be the equivalents of 6^ 7^ and 
7 7^ in Poiseuille flows, but previous work on rotat- 
ing Rayleigh-Benard convection, which can be viewed as 
analogous to our 0^0 experiments, has also shown in- 



teresting "chaotic" dynamics [22|. Last, we may wonder 
whether the forcing dependence is present in natural sys- 
tems where such reversals are observed: northern hemi- 
sphere winds zonal and blocked [l^ patterns, whereas 
Kuroshio currents might deviate abruptly from their 
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initial course: both systems exhibit complex dynamics 
with jumps between quasi-steady states. 
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